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Abstract — A composite material is prepared by joining two 
or more materials of different properties. The joined materials 
work together and give a new material with unique properties. 
Use of composite is provoked by low weight- to- stiffness and 
weight- to- strength ratios. Complex damage behaviour is 
shown by composites due to their anisotropic nature and 
heterogeneity. Thus the detailed analysis of composite 
structures is a formidable task. A joint is a structural 
connection between two or more members intended for load 
transfer. Most structures contain one or more joints. All 
structures contain joints. Joints are one of the greatest 
challenges in the design composite structures because of their 
anisotropic nature and heterogeneity, introduce high local 
stress concentrations. Damage initiation and propagation is the 
greatest concern in understanding the behaviour of bolted 
connections in composites. To support laboratory tests, a finite 
element modelling can be done to support joint design and 
predict propagation of damage. In this present study the 
analysis of a double lap joints are done using continuum shell 
elements and a progressive failure analysis was done using 
Tsai-Hill failure criteria and material stiffness degradation 
mechanism. . Progressive failure analysis was also done to 
determine the mode of failure and showed good correlation with 
the stress results. Primarily two modes of failure observed i.e. 
fibre failure and matrix failure. 

Index Terms — Composite, anisotropic, progressive failure 
analysis, mode of failure. 

I. Introduction 

A composite material is prepared by joining two or more 
materials of different properties. Use of composite is 
provoked by high specific stiffness and high specific strength 
[1]. Improved weight savings, increased fuel efficiency, 
enhanced durability, and superior structural proper-ties make 
composite materials ideal for aerospace applications [2]. 
From the library of elements available composites can be 
modeled using shell elements, continuum shell elements and 
solid elements [3]. 

A joint is a structural connection between two or more 
members intended for load transfer. All structures contain 
joints. Joints are one of the greatest challenges in the design 
composite structures because of their anisotropic nature and 
heterogeneity, introduce high local stress concentrations. 
Thus overall structural capacity is determined by the joints. 
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Several researchers have done studies on the strength of 
single lap composite bolted joints [4-6]. Effect of bolt-hole 
clearance was investigated on single-lap, single-bolt 
composite joints. Increasing clearance was found to result in 
reduced joint stiffness and increased ultimate strain in all 
tested configurations. In single-lap joints, clearance caused 
three-dimensional variations in the stress distribution in the 
laminate. These variations dependent on the lay-up sequence. 
A highly efficient user-defined finite element model and 
empirical expressions were developed to determine the 
bolt-load distribution in large-scale composite structures [7, 
8 ]. 

Damage initiation and propagation is the greatest concern 
in understanding the behavior of bolted connections in 
composites. To support laboratory tests, a finite element 
modeling can be done to support joint design and predict 
propagation of damage. Failure modes and trends in material 
response evaluated to assess the progression of failure in 
composite joints. Various progressive damage mechanisms 
are a) continuum damage mechanics (CDM) or material 
properties/stiffness degradation method (MPDM) all forms 
of damage is represented as local stiffness reduction in 
individual elements. Poisson’s ratios are not degraded and 
only the Young’s moduli and shear modulus are modified for 
a failed element, b) Discrete damage modeling (DDM) in 
which matrix cracks and delamination are explicitly 
introduced into model as displacement discontinuities, which 
they create, c) X-FEM formulations, degrees of freedom are 
added to elements along the crack surface to describe the 
displacement discontinuity, d) Cohesive elements or the 
element failure method (EFM) model formulation was used 
for crack opening [9, 10]. Prediction of the failure carried out 
using various failure criteria such as Hashin S, Tsai-Hill and 
Tsai-Wu failure theory. The results obtained were compared 
and plotted against some available experimental findings 
[11,13]. 

In this present study, validation procedure was carried out to 
determine the accuracy of SC8R continuum shell elements 
and to verify the modeling strategy. The analysis of a double 
lap joints were carried out using continuum shell elements 
and a progressive failure analysis was done using Tsai-Hill 
failure criteria and material stiffness degradation mechanism. 

II. Failure Criteria 

A successful design requires efficient and safe use of 
materials. Composite materials have many mechanical 
characteristics that are different from those of more 
conventional engineering materials. Composite materials are 
inhomogeneous (i.e. constitute non-uniform properties over 
the body) and non-isotropic (orthotropic or more generally 
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anisotropic). An orthotropic body has material properties that 
are different in three mutually perpendicular directions. Have 
three mutually perpendicular planes of material property 
symmetry. Thus the properties depend on orientation at a 
point in the body. Isotropic materials mainly have two 
strength parameters such as normal strength and shear 
strength. Failure is initiated for an isotropic material if any of 
the parameters is greater than the corresponding ultimate 
strengths. 

Theories were developed to compare the state of stress in a 
material to failure criteria. The two failure theories used are 
Tsai-Hill Failure Theory and Tsai-Wu Failure Theory, in 
which the strength parameters ( Xt, Yt, Xc, Yc and S) are 
determined through experiments and stress induced (SI 1, 
S22 and SI 2) are results obtained from Finite element (FE) 
model. 

Xt -Tensile strength in X direction 
Yt - Tensile strength in Y direction 
Xc- Compressive strength in X direction 
Yc- Compressive strength in Y direction 
S- Shear strength 

S 1 1- Stress induced in principal direction 
S22-Stress induced in transverse direction 
SI 2-Shear stress induced 

a) Tsai-Hill Failure Theory 


h 


n2 CC o 2 r< 2 

^11 211^22 . ^22 *12 

X 2 X 2 Y 2 S 2 


( 1 ) 


If I F >1, failure have occurred. 





L=100 mm 


Fig. 1. Configuration 


A. Material Properties 

The materials used for the composite laminate are carbon 
-epoxy, Glass- epoxy and steel. The material properties are 
shown in the Table 1 and Table 2. 

Table 1. Normalized material properties of carbon -epoxy 
and Glass epoxy 


Property 

Carbon-epoxy 

Glass-epoxy 

E l /E t 

16.63 

2.47 

VlT 

0.31 

.229 

Glt/Et 

0.67 

0.25 

x L /x T 

2.03 

1.70 

Y l /Xt 

0.04 

0.3 

y t /x t 

0.09 

0.45 


Table 2. Material properties of steel 


Property 

Steel 

E l (MPa) 

200000 

VLT 

.31 


b) Tsai-Wu Failure Theory 

h = ^1^11 + F 2 S 22 + F n S h + F 22 S 22 + F 66 S 12 + ^ F 12 S U S 22 
( 2 ) 

If Ip >1, failure have occurred. 


B. Composite layup configuration 

The total thickness of the composite layup is 4mm and the 
layup sequence is as shown in the Table 3. The composite 
plate is symmetric about mid layer thus only half thickness is 
been considered for analysis. Continuum shell elements are 
used to mesh a composite layup. Figure 2 shows the ply stack 
diagram of composite plate. Figure 3 shows the orientation of 
ply with respect to loading direction. 

Layup Sequence is: [90/0/-45 /45/0/90/0- G 1/2 ]s 


Where, 


L + _L 

x, x c 


F n =- 


1 


*,X C 


F 66~ — 

s 2 



F 2 2 - 


1 

Wc 


F n~f yiaa 

*/ - constant default 
value is zero 


III. Model Configuration 

In this configuration the width of the composite plate is 
changed and the effect of change in width is been 
investigated. The plate is subjected to a load of 14 kN. The 
centre to centre distance between bolts is taken as 100mm 
and the edge distance is 15mm as shown in Figure 1. The 
diameter of bolt is 10mm. The widths of the composite plate 
are 50mm, 40mm, 35mm, 20mm and 15 mm. 


Table 3. Layup configuration of composite plate 


Layer 

No. 

Material 

Thickness 

(mm) 

Fiber 

orientation 

(degree) 

1 

Carbon- epoxy 

0.3 

90 

2 

Carbon- epoxy 

0.3 

0 

3 

Carbon- epoxy 

0.3 

-45 

4 

Carbon- epoxy 

0.3 

45 

5 

Carbon- epoxy 

0.3 

0 

6 

Carbon- epoxy 

0.3 

90 

7 

Glass - epoxy 

0.3 

0 

8 

Carbon- epoxy 

0.3 

90 

9 

Carbon- epoxy 

0.3 

0 

10 

Carbon- epoxy 

0.3 

45 

11 

Carbon- epoxy 

0.3 

-45 

12 

Carbon- epoxy 

0.3 

0 

13 

Carbon- epoxy 

0.3 

90 
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Fig. 2. Ply stack plot of composite plate. 



Fig. 3. Orientation of ply with respect to loading direction 


considered. Geometric and material nonlinearity were 
included in the model. The third and final step in the 
progressive failure analysis is to apply the material 
degradation model to the failed material points. The material 
properties are degraded based upon the damage mode. The 
progressive failure analysis is implemented in Abaqus. The 
process is invoked at each material point of an element to 
evaluate the failure criterion. When failure is detected, the 
degradation model is applied accordingly. In this model, the 
material stiffness E lls E 2 2 and G i2 are instantaneously 
reduced by 1000. 


V. Modeling 

The FEM model is as shown in the Figure 5. The model 
has symmetry about X, Y and Z direction, thus only l/8th of 
the configuration is analysed. The Figure 6 shows the meshed 
model configuration used for analysis in abaqus software. 


IV. Progressive Failure Analysis 

Failure of composite structures is a progressive series of 
events. It often starts as a tiny crack between the fibres and 
matrix. These cracks reduce the stiffness of the composite. 
Capturing stress redistribution is the key to realistic 
simulation of failure in composite structures. Progressive 
failure analysis is done on the same configuration at failure 
load. Progressive failure analysis helps us in knowing the 
mode of failure. The damage in composite structures is 
generally a combination of matrix cracking, fibre breakage in 
tension and compression, and delamination. The first two 
damage modes, matrix cracking and fibre breakage in tension 
and compression are considered. 



Fig. 4. Flow chart showing process of progressive failure 
analysis 

The procedures for predicting the growth of the damage 
path are developed using the progressive failure analysis 
methodology implemented within finite element analysis. 
The progressive failure analysis methodology generally 
consists of three steps Figure 4 shows the flow chart of 
methodology of progressive failure analysis: a) calculating 
the lamina stress { Stresses computed in principal, transverse 
and shear directions} b) Estimating failure index and c) 
degrading the material stiffness in the failed elements to 
represent damage. In the study intra-laminar failure modes 



r*m 

t?T 


Tff 


hi 

IT “Ti 


TTT 



Fig. 5. Configuration model in abaqus 



Fig. 6. Configuration meshed model in abaqus 
VI. Results 


A. Validation 

The accuracy of any FE model is dependent on the accuracy 
of the geometry, the type and number of elements used, and 
the material property model. Validation is done in this study 
to check whether the SC8R elements used produce required 
results for the composite layup and to check the modeling 
strategy. For that purpose a problem done by Buket Okutan 
[14] is selected and FE modeling was carried out. The results 
are verified with the results obtained from previous study. 

Geometry: A rectangular composite plate has length L, 
thickness t and width W with a hole of diameter D. The hole 
is at a distance E from the free edge of the plate . The 
configuration of composite plate is shown in Figure 7. In the 
study [14] it was observed that for the [0/90/0] s laminate, 
failure modes were found as bearing mode when the E/D 
ratio is greater than 3. Thus E/D ratio 4 was taken for 
validation purpose. 
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E=4flmm L=8Gmm 
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Fig. 7. Geometry of specimen 

Material: Table 4 shows the material properties of 
glass-fiber/epoxy composite. Stacking sequence of 
composite plate is [0/90/0] s 

Load: A tensile load is applied at the hole free edge of the 
plate resisted by the pin. 

FE model and boundary condition: Thus composite plate 
modeled using abaqus software gives required results. 
Composite plate was modeled using continuum shell element 
(SC8R). Widths modeled are 20mm, 30mm, 40mm and 
50mm. Figure 8. shows FE model done for the joint 
configuration. l/4 th of configuration was modeled due to 
symmetry in Y and Z direction. Symmetric boundary 
conditions was been applied. Tensile load is applied at the 
hole free end of composite. In the problem pin was assumed 
to be rigid and thus not modelled. The degrees of freedom 
were arrested in the quarter portion of bolt hole to simulate 
support conditions. 

Table 4 Properties of glass -fiber/epoxy composite 


Longitudinal modulus El (MPa) 

44,000 (+560) 

Transverse modulus E2 (MPa) 

10,500 (±420) 

Shear modulus G12 (MPa) 

388045 (±360) 

Poisson’s ratio v i2 

0.36 

Longitudinal tension Xt (MPa) 

800 (±59) 

Longitudinal compression Xc (MPa) 

350 (±42) 

Transverse tension Yt (MPa) 

50 (±4.35) 

Transverse compression Yc (MPa) 

125 (±9.34) 

Shear strength S (MPa) 

120 (±15.28) 
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Fig. 8. FE model done for the joint configuration 


Results: The obtained results (Figure 9) showed the 
variation of bearing strength with w/d ratio The results show 
good correlation and thus the results obtained are validated. 
Thus SC8R elements can be used to model the composite 
layup. In the present study bolts are also modeled to replicate 
the contact property in the real problem. 



Fig. 9. Variation of bearing strength with w/d ratio 

The FEM model is as shown in the Figure 5. The model 
has symmetry about X, Y and Z direction, thus only l/8th of 
the configuration is analysed. The Figure 6 shows the meshed 
model configuration used for analysis in abaqus software. 


B. Stress Results 


Figure 10 shows Sll stress distribution in each layer of 
composite layup of single pin configuration of width 15mm. 
It can be seen from Figure 10 that, stress is concentrated at 
the bolt hole. Thus the circumference of bolt hole is set as the 
region of interest. The Figure 11 -Figure 13 shows the 
variation of stress, i.e. Sll, S22 and S12 around the bolt hole 
in each layer of the composite plate. The Sll, S22, S12 
stresses are evaluated and are used for the calculation of 
failure index. Tsai Hill (TSAIH) and Tsai Wu (TSAIW) 
failure criterion is used to determine the failure index of each 
layer and the layer in which first failure (first ply failure) 
occurred is determined. The component that is responsible 
and the failure load are also determined. Table 5 & Table 6 
shows the failure assessment details using Tsai-Hill and 
Tsai-Wu failure theories. 
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Fig. 10. Stress Sll in each layer of 15mm configuration 


L n y up s e quenc e (90 0 -45 45 0 90 0 [Ol 2] ) sym 
3.00E+03 
§2.00E+03 
g 1.00E+03 

i/T 

lo.OOE+OO 

VI 

-1.00E+03 
-2.00E+03 

Circumferential angle (Degree) 

Fig. 1 1 . Variation of stress Sll around the bolt hole in each 
layer of the composite plate for single pin configuration 




Fig. 12. Variation of stress S22 around the bolt hole in each 
layer of the composite plate for single pin configuration 



Fig. 13. Variation of stress S12 around the bolt hole in each 
layer of the composite plate for single pin configuration 


Table 5. Failure assessment details of configuration by Tsai-Hill failure criteria 


Width 

(mm) 

Load 

(kN) 

Tsai Hill 

failure 

index 

Initial 

Failure 

layer 

Number 

Initial Failure 
layer 

orientation 

Location of failure 
with respect to 
circumference angle 
(degree) 

Failure 

component 

15 

3.5 

3.286 

6 

90 

95.58 

S22 

20 

3.5 

2.161 

6 

90 

95.58 

S22 

35 

3.5 

1.740 

4 

45 

52.23 

Sll 

40 

3.5 

1.726 

4 

45 

52.23 

Sll 

50 

3.5 

1.72 

4 

45 

52.23 

Sll 


Table 6. Failure assessment details of configuration by Tsai-Wu failure criteria 


Width 

(mm) 

Load 

(kN) 

Tsai Wu 
failure 
index 

Initial 

Failure 

layer 

Number 

Initial Failure 
layer 

orientation 

Location of failure 
with respect to 
circumference angle 
(degree) 

Failure 

component 

15 

3.5 

3.421 

6 

90 

95.58 

S22-square 

20 

3.5 

2.270 

6 

90 

95.58 

S22-square 

35 

3.5 

1.991 

4 

45 

52.23 

SI 1-square 

40 

3.5 

1.972 

4 

45 

52.23 

SI 1-square 

50 

3.5 

1.971 

4 

45 

52.23 

SI 1-square 
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C. Progressive Failure Analysis 
a) 15 mm width 

Table 7. Progressive failure analysis of configuration 15mm 


No 


Failure 

load 

(kN) 


Layer 

failing 

Sequence 


Figure 


1 


2 


3 


4 


5 


1.095 


1.095 


1.095 


L6,L1 


L6,L1,L3 


L6,L1,L3, 


L4 


1.095 


L6,L1,L3,L4 

,L5,L2 


1.095 


L6,L1,L3,L4 

,L5,L2,L7 



6 


7 


8 


1.095 


1.095 


1.095 


L6,L1,L3,L4 

,L5,L2,L7 


L6,L1,L3,L4 

,L5,L2,L7 


L6,L1,L3,L4 

,L5,L2,L7 



Failure of laminate is assumed to occur when the element 
degradation reached up to the plate edge. Progressive failure 
analysis results of 20mm configuration are similar to Table 7. 
b) 35mm width 


Table 8 Progressive failure analysis of configuration 35 mm 


No 

Failure 

load 

(kN) 

Layer 

failing 

Sequence 

Figure 

1 

2.03 

L4 







2 

2.03 

L4,L1,L2,L5 

,L6 






3 

2.03 

L4,L1,L2,L5 
,L6, L3,L7 




4 

2.03 

L4,L1,L2,L5 
,L6, L3,L7 





5 

2.03 

L4,L1,L2,L5 
,L6, L3,L7 






Failure of laminate is assumed to occur when the element 

degradation reached up to maximum displacement limit. 

Similar failure results Table 8. are obtained for 40mm and 

50mm configurations. 

VII. Inferences 

• The stresses are concentrated at bolt-hole regions. 

• The location of maximum stress concentration depends on 
the fibre orientation of each layer. 

• The stress plot around the circumference of hole follows a 
particular trend for a fibre orientation. 

• The stress concentration is larger for carbon epoxy layer 
than glass epoxy layer. 


Table 9 describes the mode of failure of composites. There 
are two modes of failure mainly fibre failure and matrix 
failure. 

Table 9 Mechanics of failure of composites 



Stress component 

Mode of failure 

l 

Sll 

Fibre failure 

2 

S22 

Matrix failure 


Table 10 shows the summary of progressive failure 
a nalysis. Figur e 14 . variation o f failure load w ith width. 



0 2 4 6 

w/d ratio 


Fig. 14. Variation of failure load and width of configuration 
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Table 10. Summary of analysis of Single pin configuration 


Width 

(mm) 

Failure 

Load 

(kN) 

Initial 

Failure 

layer 

No: 

Initial 

Failure layer 
orientation 

Location of failure 
with respect to 
circumference 
angle (degree) 

Failure type 

REMARK 

1.5 

1.095 

6 

90 

95.58 

Tensile Failure 

S22 stress is Tensile 

2.0 

1.62 

6 

90 

95.58 

Tensile Failure 

S22 stress is Tensile 

3.5 

2.03 

4 

45 

52.23 

Bearing Failure 

S 1 1 stress is compressive 

4.0 

2.04 

4 

45 

52.23 

Bearing Failure 

SI 1 stress is compressive 

5.0 

2.045 

4 

45 

52.23 

Bearing Failure 

SI 1 stress is compressive 


VIII. Conclusion 

The use of composites in load bearing structures is 
primarily motivated by high specific stiffness and high 
specific strength. The stress distribution depends on the 
layup sequence and materials used. First ply failure occurs 
when the first ply or ply group fails in a multidirectional 
laminate. Progressive failure analysis was carried out to 
determine the mode of failure and showed good correlation 
with the stress results. 

• The stresses are concentrated at bolt-hole regions. 

• The stress concentration is larger for carbon epoxy layer 
than glass epoxy layer. 

• Modes of failure considered are fibre failure and matrix 
failure. 

• When, w/d < 2 : Failure type is tensile failure 

• When, w/d>2 : Failure type is Bearing failure 

• Failure load increases as the width of plate increased. 
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